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Abstract 

We consider an inverse boundary value problem for the hyperbolic 
partial differential equation 

n 

{-id t +A (t,x)) 2 u(t,x)-Y^ (-id x .+A j (t,x)) 2 u(t,x)+V(t,x)u(t,x) = 

i=i 

with time dependent vector and scalar potentials (A = (Aq, . . . ,A m ) 
and V(t, x) respectively) on a bounded, smooth cylindric domain 
(— oo, oo) x f2. Using a geometric optics construction we show that the 
boundary data allows us to recover integrals of the potentials along 
'light rays' and we then establish the uniqueness of these potentials 
modulo a gauge transform. Also, a logarithmic stability estimate is 
obtained and the presence of obstacles inside the domain is studied. 
In this case, it is shown that under some geometric restrictions similar 
uniqueness results hold. 

1 Introduction 

Let Q be a bounded domain in M n , n > 2, consider the hyperbolic equation 
with time dependent coefficients 

n 

(-id t + A (t,x)) 2 u-^~2(-id Xj +A j (t,x)) 2 u + V(t,x)u = inixfi, (1) 
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where V(t,x), Aj(t,x), < j < n, are smooth functions vanishing when 
{\x\ > R} for some R > 0. The smooth vector field A(t, x) = (A (t, x), . . . , A n (t, x)) 
is called the vector potential, the function V(t, x) is called the scalar potential 
and equation ([1]) is often referred to as the relativistic Schrddinger equation 
(see [21]). 

For the above differential equation we impose the initial and boundary 
conditions 

u (t, x) = d t u(t, x) = for t « (2) 

u{t,x) = f (t,x) onRxffl, (3) 

where / is a compactly supported smooth function on R x dQ. Solutions 
to ([1]) satisfying ([2]) and ([3]) are unique and we can define the Dirichlet to 
Neumann operator by 

A(/) : = (d v + iA(t,x) ■ v)u{t,x) (4) 

where u is the solution of ([I])-©, v is the exterior unit normal to dQ and 
we have set A(t,x) = (Ai(t, x), . . . , A n (t, x)). The Inverse Boundary Value 
Problem is the recovery of A(t,x) and V(t, x) knowing A(/) for all / G 
Q° (R x 80) . 

Inverse problems is a topic in mathematics that has been growing in 
interest for the past decades, in part, due to its wide range of applications, 
from medicine to acoustics to electromagnetism just to mention a few (see 
for instance p3J for some of the latest tools and techniques employed in the 
solutions of these problems). In the case of the hyperbolic inverse boundary 
value problem (l)-(4) with time independent coefficients, a powerful tool 
called the boundary control method, or BC-method for short, was discovered 
by Belishev (see |3J). It was later developed by Belishev, Kurylev, Lassas, 
and others Ql7|.[18]). and more recently a new approach to this problem 
based on the BC-method was developed by Eskin in ([5], [6]). On a similar 
note, Stefanov and Uhlmann established uniqueness and stability results for 
the wave equation in anisotropic media (see [22] and [30] for a survey of these 
results). 

Nevertheless, the case of time dependent coefficients has seen very little 
progress in recent years. In the case of the vector potential being identically 
equal to zero (^4 = in (1)), Stefanov [25] and Ramm-Sjostrand [23], have 
shown that the Dirichlet to Neumann map completely determines the scalar 
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potentials. More recently, Eskin [7] considered the case of time-dependent 
potentials that are analytic in time. The analiticity of the time variable is 
related to the use of a unique continuation theorem established by Tataru in 
[28] . In this paper we eliminate the restriction on the analiticity in the time 
variable and not only we extend the uniqueness results in [25], [23J, but we 
also establish a logarithmic stability estimate for the case when the vector 
potentials are compactly supported in space and time. 

We also study the problem with obstacles inside the domain and show 
that under some geometric considerations similar uniqueness results hold. 
The presence of these obstacles in the domain may lead to the Aharonov- 
Bohm effect. This problem was considered by Nicoleau and Weder in the 
context of the inverse scattering (see [20] and [21] respectively), and by Eskin 
[8] in the context of the inverse boundary value problem for the Scrhodinger 
equation. 

This work is structured as follows. In section [2] we introduce the notion 
of gauge equivalent for a pair of vector and scalar potentials and we make 
some remarks about uniqueness. In section [3] we construct geometric optic 
solutions (GO for short) for equation ([T]) satisfying the set of initial con- 
ditions ([2]). In section H] we establish a Green's formula for these types of 
problems and show that the light ray transforms of gauge equivalent poten- 
tials agree. In section [6] we prove uniqueness of the potentials in the case 
where no obstacles are allowed inside the domain Q. In section [7J based on 
the works of Isakov [T3], Isakov and Sun [T3] and more recently Begmatov 
[2], we establish a stability result of logrithmic type for the particular case 
when the potentials are compactly supported in both space and time. Finally 
in section [8] we consider the problem when one or more convex bodies are 
allowed inside the domain (by imposing a geometric restriction on the layout 
of these obstacles). 

2 Gauge equivalence 

The ultimate goal in most inverse boundary value problems is the recovery 
of the coefficients of a partial differential equation, however for application 
purposes this recovery is meaningless unless it can be done in some sort of 
'unique' way. In our case this type of uniqueness is obtained modulo a gauge 
transform. 
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Definition 2.1. We say that the vector and scalar potentials (A(t, x), V(t, x)) 
and (A(t, x), V'(t, x)) are gauge equivalent if there exists g(t,x) G x 
Q)) such that g(t, x) ^ on M x Vt), g = 1 on R x <9f2 and 

A'(t,x) =A(t,x) - —^—Vt. x g{t,x) 
g(t,x) ' 

V(t,x) =V(t,x), 

where V t , x '■= ipu d x ) = (d t , d xi , . . . , d Xn ) is the (n+ l)-dimensional gradient. 
The mapping (A, V) — > (A, V) is called a gauge transform. 

The definition above includes the more general case when obstacles are 
present inside the domain. When Q is simply connected (no obstacles), the 
gauge g has the particular form g(t, x) = e % ^ 1 ^ where cp(t, x) E C°°(R x fl). 
Then — V( tjX )g(t, x) = V(t >a; )</?(t, x) and we see that two vector potentials 
are gauge equivalent if their difference is the gradient of a smooth function. 
The following proposition tells us that recovery of the potentials can only be 
done up to a gauge transform. 

Proposition 2.1. Ifu(t,x) is a solution of (QP-([3j] and g(t,x) is as in defi- 
nition \2. then v(t,x) = g(t,x)u(t,x) satisfies 

n 

(-id t + A'o(t,x)) 2 v - + A'j{t,x)) 2 v + V'(t,x)v = in R x Q 

(5) 

v = d t v = for t < Ti 
v = /<?L an onR x dfl. 

with {A ', V) and {A, V) gauge equivalent. 

In addition if A' is the Dirichlet to Neumann operator associated to (TJP, then 

A'(y |Rxdtt) = <?|]RxdnA(u|]Rxan) (6) 
i.e., A' = A since g\nxdn = 1- 
Proof. Setting xq = t we see that for < j < n, 

(-id Xj +A' j (t,x))v(t,x) =g(t,x) \ -id x . + A'^t, x) - Z d Xj g(t, x)\ u(t, x) 
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Choosing A'-{t, x) so that A'- = Aj + ~d Xj g for < j < n, we get 

(-id x . +A' j (t,x)) 2 v(t,x) = (-id^+Afax)) (g(t,x) (-id x . + Aj(t, x)) u(t, x)) 



thus 



g(t,x) (—id Xj + Aj(t,x)) u(t,x), 



(-id t + A' (t,x)) 2 v i9 *, + A'^xtfv + V'(t, 



XjV 



n 

g(t, x)((-id t + A (t, x)) 2 u-J2i- id *i + Ml x )f u + v & x ) u ) = 

3=1 

as u is a solution of ([1]). Also notice that since g is smooth and u satisfies 
(TJD and © we have for t « 

v(t, x) = u(t, x)g(t, x) — 
d t v(t, x) = u(t, x)d t g(t, x) + d t u(t, x)g(t, x) = 0, 



similarly v(t, x) 
simply notice that 



ixdn 



(g(t,x)u(t,x))\ 



ixdn 



fg 



ixdW 



To conclude we 



A' v 



ixan 



(d„(gu) + %A' ■ is(gu)) 



ixdn 



(d u g) u + g (d u u) + i(A + ig l d x g) ■ u(gu) 



g(d u u + i(A ■ v)u) + {(d v g)u - (d u g)u) 



ixdn 



Rxdn 



ixdn 



ixdn 



An 



ixdn 



□ 



If the above equality holds we shall say that the Dirichlet to Neumann 
maps A and A' are gauge equivalent. Summarizing, we have shown that if 
the vector and scalar potentials are gauge equivalent then the Dirichlet to 
Neumann maps are equal. In the following pages we shall attempt to prove 
the converse, roughly speaking: If for a pair of vector and scalar potentials 
the Dirichlet to Neumann operators associated to the hyperbolic equation (QP- 
(T3|) are equal, then so are the vector and scalar potentials. 
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3 Geometric optics 

For the hyperbolic problem (JT])-([3]) we shall attempt to construct geometric 
optics solutions supported near light rays. In order for us to achieve this goal 
we consider solutions having the form 

u(t, x) = e »(^-) Vj §mr + v(N+1) ^ *)> «> e S"~\ k e R. (7) 

For u as above we have 

JV 

(-id t + A )u = e ik ^^ (-id t + Ao + ik)(Y] 7T7TV- + e-^—V^ 1 )) 

p=0 v ' 

(8) 

applying the above identity twice to a solution of (JT]) we get 
(-id t + A )\ = (-id t + A ) (e*(*—) ( - id t + A + ik) (J2 J^y p + 



e -ik(t-wx) v (N+l) 



» 



JV 



e -ik(t-wx) v (N+l) 



this is 

( - id t + A ) 2 u = e ih ^ x H( - id t + A ) 2 + 



JV 

2ik( - id, + - *») (£ JL- + e-*«— V"+'> 



p=0 

Since a similar formula holds for 1 < j < n, we obtain 
( - id Xj + A 3 )\ = e ik ^ x ) (( - + ^) 2 + 

JV 



Sto, ( - + A,) - *;S 2 ) (E(^)? + ^ 



thus equation becomes 

n 

= Lu= e ik ^ x ) (( - idt + A ) 2 - ( " i9 Xj + A,) 2 + 

n 

+ 2ike lk ^^ (( - id* + A ) + ( - z^, + 

n 

i=i 

= Lu = e ik ^-^ (L + 2ik£)v, (9) 



\2' 



where we have set 

JV 

' (2ik)P 



= E^r + ^^^H^) (10) 

p=0 



L = ( - id t + A (t, x)) 2 -Ys(- i9 Xj + Aj{t, x)) 2 + V(t, x) (11) 



n 



£ = (-i^ + A (t,x)) +^2^j{-id Xj +Aj(t,x)). (12) 

Plugging in the expression for v into (J9j) we obtain 

= (2ikC + L) + ^ , + - + ^„ + e -C—M™>) , (13) 
which in turn can be rewritten as 

(2ik)£v + (Cv! + Lv ) + -— (£i; 2 + LuJ + • ■ ■ + 

(14) 

where we have used the identity (2ikC+L) ( e -^(*-<^y^+i)) = e -ifc(t-w*)_^,(JV+i)_ 

Notice that a solution of ([I])-® can be found by solving the iV + 1 trans- 
port equations 

Cv = 0, Cvj = -Lvj-x, l<j<N (15) 
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with initial conditions supported near a neighborhood of the light ray 7 = 
{(t',x') + s(l,u)) : (t',x') J_ (l,o;), s G R} (we assume that 7 intersects 
the plane t = T\ outside of the cylinder R x O; as well as the second order 
equation 

(2i£;) J 

with initial and boundary conditions 



Lv {N+1) = -7K^Lv N (16) 



v {N+1 \t,x) = for* = Ti 

d t v (Ar+1) (t,a;) = for t = T x 

v (N+1) (t, x) = for t >Ti, x £ dtt. 

The above differential equation has a unique solution. Moreover if h is 
the right hand side of (ITS]) we have (see for instance Isakov [TJ], pp. 185) 
that if Tx < t < T and jfe > 1 

\\v {N+1 \t, )|| i2(n) <C||/i|| i2((Tl)T)xn ) 

Thus, we have shown that we can find a solution u = e ik( - t+UJ - x \v + Oik- 1 )) 
of ([1]) satisfying the set of initial conditions ([2]). Let us now examine the first 
term in fTTUT) by solving the transport equation 

n n 

= Cv (t,x) = s ^u j d Xj v (t,x) + i^2ujAj(t, x)v (t, x) (18) 

j=0 j=0 

where we have set o;o = 1 and d Xo = d t . 

Equation (Tl8l) is a first order transport equation that can be solved by 
the method of the characteristics or by performing a change of variables that 
turns the PDE into an ordinary differential equation. Either way, the solution 
we obtain is given by 

/ i (t+bJ . x) n ^ ^ \ 

t>o(t, x) — x{t, x) exp I — i I 2_, u jAj{t' + s, x + su) ds, J (19) 

V j=o J 

where (t',x') = (t,x) — |(t + 0; ■ x)(l,u) is the projection of (t,x) into II( ljtJ ), 
the hyperplane perpendicular to (1, 0;) and \ is any function that is constant 
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along the direction given by (1, u) and whose support is contained in a neigh- 
borhood of the light ray 7 = {(£', x') + s(l, u>) | s G M} (in general x can be 
complex valued but for our purposes we will assume it is real valued). 

Summarizing, we have been able to construct a solution of ([I])-© having 
the form: 

u(t, x) = e *(t-«-s)-«ii (*,*;«) ( x (f, x >) + Oik' 1 )) , (20) 

where 

Rx(t,x;u)= / ^ WjAjtf + s,x' + soj) ds (21) 

In a similar way one can obtain geometric optics solutions for the back- 
wards hyperbolic problem 

Lv = in(-oo,T 2 )xft 
v = d t v = for t = T 2 

In the following section we will derive a Green's Formula for these kinds 
of hyperbolic operators and will use the Geometric Optics representations to 
conclude that the Dirichlet to Neumann data determines the vectorial and 
scalar ray transforms of the potentials along 'light rays' (this is, rays that 
make a 45 degree angle with the hyperplane t — 0). 



4 Green's formula 

This technique has had a lot of success in the context of inverse problems, in 
particular for the case of elliptic problems, the fundamental paper of Sylvester 
and Uhlmann [27] has been a source of inspiration for several other unique- 
ness results (see also Isakov's review paper [13] for more information on this 
subject). 

For T\ and T 2 two real numbers with T\ < T 2 we consider the forward 
and backward hyperbolic equations 



u 



L\U = 
= d t u = 
u = f 



in [Tx,T 2 ] x n 

for t = T x 

on [T\,T 2 ] x Oil 



L* 2 v = 
= d t v = 
v = g 



in [Ti,T 2 ] x ft 

for t = T 2 

on [T\,T 2 ] x <9ft, 
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where we have set 

n 

L\ = L{A^\ VV>) = {-id t + A$\t,x)) 2 id *i + A f^ x )f + V {1) (t,x) 



3=1 



L* 2 = L(A&>, W) = {-id t + Ag\t,x)) 2 id *i + A f^ x )Y + V^(t,x) 

3=1 

and let us assume that the Dirichlet to Neumann operators 

Ai(/) = (d v + iu ■ A^(t,x))u(t,x)\ [TuT2]xdn (22) 
M9) = {d„ + ivAW(t,x))v(t,x)\ [TuT2]x9n (23) 

equal on (Ti,T 2 ) x dQ, i.e., Aif = A 2 f for all / smooth and supported on 
the set (Ti,T 2 ) x dQ. 

Remark: Notice that for the operator L\ we associate the Dirichlet to 
Neumann map 



A* 2 (g) = {d u + tvA(i)(t,x))v(t,x)\ Rxm 

and that our main assumption is A\ = A 2 on R x dQ. This is no mistake as 
later on we will show that our notation is justified as the L 2 adjoint of A 2 is 
indeed A 2 . 

Denoting by ( , )[ Ti ,t 2 ]xQ, ( , )n and ( , )[ Tl ,T 2 ]xdn the L 2 inner products 
in [Ti,T 2 ] x Q, Q and [Ti,T 2 ] x <9f2 respectively we obtain the following 
integration by parts formulas for A^ 



((-id t + A^)\v) [TuT2]xn = <( - id t + A^)u, (-id t + A^)v) [TuT2]xn 



i((-id t + A^)u(t,-),v(t,-)) ( 



T-2 



where v = (u^\...,u^) is the exterior unit normal to <9f2 and u, v are 
solutions of the forward and backward hyperbolic equations respectively. 
In view of the initial conditions we obtain 



((-^t + A ( i ) ) 2 u,v) [TuT2]xn = ((-id t +A^)u, (-^ t + A ( i ) )v) [TuT2]xn . (24) 
Also for j = 1, . . . , n we have for A^ 



<( " + Af)\v) [TuT2]xn = <( - id Xj + A?)u, ( - id Xj + A^)v) [TuT2]xn 

-i((-id Xj +A?)u^\v) [TuT2]xm , 
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hence 

n 



[Ti,T 2 ]xn 



E((-^ + 4 1 V(-^ + 4>) [Ti 



[Ti,T 2 ]xO 



(25) 



where / = u\ [TuT2]xdn and g = v\ [TuT2]x9n . 
Similarly for LjL we have 



(u, (-id t + AV) 2 v) [TuT2]xn = ((-id t +AP)u, (-id t +A®)v) [TiT2]xn (26) 
and 



,T 2 ]xQ 



£ <( - + 4 2 >> ( - ^ + A r>>[T 11 T 2] xn 

(27) 

Combining expressions (I24p - (I27I) and recalling that it and i> are solutions 
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of the forward and backward hyperbolic problem respectively we obtain 

= (L lU ,v) [TuT2]xn - (u,L* 2 v) [TuT2]xn 



<(-^ + 4V(-^ + 4» [Tl , r2]x n 



h 



-((-id t + A^)u,(-id t + A^)v) 



[Ti,T 2 }xQ 



~ E <( " id *> + ( - ^ + 4>W 2 ]xn 

i=i 



h 



+ E (( - ^ + ( - ^ + ^Mp^n 



-v 

/4 



+ (^. ? ') [ T 1 ,T 2] x,-( 1/(2) «^) [ T ll T 2 ,x, 

+ </. A 2(^)> [ T 1 ,T 2] xan - (HfM^ ]xga (28) 

Let us now study the terms Ji, I2, h and Z4 appearing in the above formula. 
For Ii and I 2 we have 



h + / 2 = < ( - idt + 4>, ( - ^ + 4r )^> [ T 1 ,T 2] xn 

-<(-^ + 4 2) )^(-^ + Af)t;> [Tiir2]xn 



- id t u, -i^) Kxn + < " ^Mp^xn + -^> [ T 1 ,T 2 ]xn 

-(^0 V _l ^ ?, )[r 1 ,T 2 ]xSi _ (^0 V ^0 M[T l5 T 2 ]xfi' 

this is 

h + h = {(A$> - Af)u, -s t v) (Ti ^ ]xn + - 4 2 ») ( - ia t «),«) 1TiTi|)<a 

+ <((4 1| ) 2 -(4 2| ) 2 )«,«V„ T2lxn (29) 
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and a similar computation shows that 
h + h = £ < (A« - Af -id Xj v) [Ti T2]xn 

3=1 

n 
3=1 

n 

+ E<((4 1, ) 2 -(lf , ) 2 )^«>,r,.r J „ n - (30) 

3=1 

Combining f l28|) . (129]) and fl30|) we obtain the Green's formula: 

( A l(/)^> [Tl ,T 2 ]x0n - (/' A 2^)>[T 1 ,T 2 ]xM = 

i=o 

n 

+ E^<((^ 2) ) 2 " (^fl^W^n " <^W]xn ( 31 ) 
i=o 

where we have set x = t, A, = Aj 2) - A^ for < j < n, V = 1/ (2) - V {1 \ 
r = — 1 and rj = 1 for 1 < j < n. 

At this point it is convenient to notice that if we take the vector and 
scalar potentials in the forward and backward hyperbolic equation to be the 
same (i.e., (A {1 \V^) = (A {2) ,V^)), then we get from (j2g]l 



( A (f),9) [Tl ,T 2 ]xdn - (fi A *(9)) [Tl ,T 2 ]xdn = 



proving that A* = d u + iv ■ A(t, x) is the L 2 adjoint of A = d v + iv ■ A(t, x). 

Since we are assuming that the Dirichlet to Neumann maps for the for- 
ward and backward hyperbolic equations agree (A 1 = A 2 on R x dQ), then 
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equation (I3TT) can be rewritten as 

= <Ai(/),^> [TiT2]xan - <A 2 (/)^) [TiT2]x ^ 

= (Mf)>9) [TuT2]xdU -(f,K(g)) [TuT2]xdn 

n 

= J2 r i{( A J u > H'V^.iyxn + <^'(-^ M )^>[T 1 ,T 2 ]xn) 

3=0 

n 

+ 5><((4 2) ) 2 - (^H-W^ - ( Vu A TuT2]xn ( 32 ) 

3=0 U 2 

where as before x = t, A 5 = Af - Af for < j < n, V = - V®, 
r = — 1 and rj = 1 for 1 < j < n. 



5 X-ray transform 

Our next step is to combine our two main tools, namely the geometric op- 
tics representation of the solutions of the forward and backward hyperbolic 
equations and the Green's formula. 

Owing to ( 120]) and (12 ip we can write geometric optics representations 
for u and v the solutions of the forward and backward hyperbolic equation 
respectively. These representations are 

u(t, x) = e W-"-x)-iRi(t,*;") x >) + o (k- 1 ) ) (33) 

^J~^ = e -ik(t-u,-x)+iR 2 (t,x;u,) x /) + Q ^-1 j j ^ (34) 

where 

R 1 (t,x;uj)= ^ujjA^it' + s,x' + su)ds (35) 

i=o 

R 2 {t,x;u)= ^oojAfit' + s,x' + su)ds. (36) 

Notice that for u as in (1331) . we have for < j < n 
d x .u = e *Kt-™)-iRi(tw) (d x . x + o (k- 1 ) + (-ikrjuj - id Xi Ri) (x + O (&T 1 )) 

= ^gifcCt-w-arJ-ifl!^) r _ ir . u . x + ©(AT 1 )) , 

14 



where we have set Uq = 1 and as before r = — 1 and rj = 1 for 1 < j < n. 
Then owing to (1341) we obtain 

{-id Xj u(t,x)) tfj~x) = - ke iMtw)-Mtw)) ( rjUJjX (t\x') 2 + Oik- 1 )) 
similarly 

u(t,x){-id Xj v{t,x)) = - k ei R rt> x ^- R ^> x ^) ( r jUjX (t', x'f + Oik^ 1 )) 
Thus, Green's formula now reads 

it x Jn j=0 

e i{Mtw)-Riitw)) dx dt + 



= Ck P [ Y^(Af\t,x)-Af{t,x)y j uj j x 2 {t\x') 



where C is a (negative) constant and "• • •" represents terms of order 0(1). 
Dividing the above expression by Ck and taking the limit as k — > +oo we 
get 

0= p j J2cuj (A?\t,x) - Af\t,xj) x^t^x^ 1 ^^-^^) dxdt. 

JT\ j_Q 

Without loss of generality (cf. Remark 3.1 in [8i]) we can assume that 
supp^ (j) C Kxfl, j = 1,2. Writing X' = (t',x') and setting A = 
(A , . . . , An) = A^ — A^ we get after the change of variables (t, x) = 
a(l,co) + X' 

/» POO n 

0= J2 u i A i ( X ' + ^ w )) ( X ') e</ "~ E ^^^ (X ' +s(1 ' w))ds da dS x >. 

Since % is an arbitrary function of X' we then conclude that 

/OO n 
OO „-_n 

da V / 

= -i ( e iJZ Z2^»i A AX'+'M)*> _ (37) 

Summarizing we have proven the following 



-» ,-=o 

"OO 

OO 
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Lemma 5.1. Suppose that the Dirichlet to Neumann operators Ai and A 2 
for the hyperbolic equations 

n 

L k u = ((-^+A« (t,x)) 2 -^ {-id Xj +A i - l \{t 1 x)) 2 +V {l) {t 1 x))u = 0, k = 1,2 

egwa/ on [Ti,T 2 ] x 90. TTien /or any light ray 

7 = {(£', x') + s(l, w) :sGl,u6 S 171-1 , (*', • (1, u) = 0}, 

£/ie vectorial ray transform of A = (A^ — A^\ . . . , A^ — An) along 7 is an 
integer multiple of lit. This is 

/oo n 
UjAj (t' + s, x' + su) ds = 2vrr (38) 
00 3=0 

for some r 6 Z. i/ere (i', x') = (i, rr) — |(i + a; • x)(l, u) and ojq = 1. 
Proof. Equation (I3TI) can be rewritten as 

e i{VA){t,x;u,) = l 

which in turn implies (I38|) . □ 

If we now incorporate the hypothesis of and A^ being compactly 
supported in x we can determine the exact value of r. This is because equa- 
tion f )38|) holds for any (£, a;; a;) G R™^ 1 x 5* m_1 and in particular, when 
t — and 1 2 1 is big enough and perpendicular to a fixed w, the light ray 
(0, x) + s(l,ou), s e K does not meet the support of ^4, hence 

/oo m 
UjAjtf + s,x' + slo) = 0. (39) 
00 j=0 

To conclude this section let us proceed to remove the condition {t! , x') ■ 
(1, to) = 0. If (t, x) is an arbitrary point in M™^ 1 , then by making the change 
of variables s = a — t+ ^ x we get 

/oo n /*oo n 

yj UjAj(t + a,x + aco) da = UjAj(t' + s, x + stu) ds, 

-00 j=0 J-00 i=0 

where (t',x') = (t,x) — (1 ; o>). Clearly this last integral equals zero by 
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6 The main theorem 



In this section we will establish several uniqueness results for vector and 
scalar potentials satisfying different growth conditions. Let us proceed first 
with the part that deals with the vector potentials in the case when the 
component functions Aj(t,x) decay exponentially in t. 

Theorem 6.1. Suppose that A(t, x) = (A (t, x), . . . , A n (t, x)) with A G C°° 
in x and t is such that for any non-negative integers a, f3 and for any < j < 
n there exist positive constants c, C a ^ such that for \t\ > t , \d^d^Aj(t,x)\ < 
Ca,/3e _c '*'. If in addition Aj(t,x) = for \x\ > R > and (3$j) holds, then 
there exists (p(t,x) G C°°(t,x) and positive constants C' a o,c' such that 

i) A (t,x) = d t ip{t,x), Aj(t,x) = d x .(p(t,x), l<j<n, and 

ii) Supp (p CRx {|x| < R}, \dfd^{t,x)\ < C; i/3 e- c 'l*l. 
Proof. By uniqueness of the Fourier transform and by fl39|) we have 




By the hypothesis on the support of Aj, 1 < j < n, we can change the 
order of integration. After writing t\ = t + s, x\ = x + su, we are lead to 

= HI e -iiti-s)T-i{*i-su>u ^ Q + cUjA^J (t h xi) dhdxxds 

j=l 

/n . 
e^+-O^A + J2^) (T,£)ds 

3=1 
n 

= S(t + u-0{a q + J2^j) (t,0, (40) 

3=1 

which tells us that the Fourier transform of Aq + Y^=i vanishes on 
11(1,0;), the hyperplane perpendicular to (l,u>). 

We claim that this Fourier transform vanishes in the complement of the 
'light cone' C = {(r, £) : \r\ > |^|} for an appropriate choice of co. To see 
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this notice that if (r, £) ^ C, then -j^j 1 has norm less than one and we can find 
u = u(t,£) G S 71 ^ 1 satisfying 

w a ™ = -® (41) 

Clearly with this choice of w we have r + u>(t, £) • £ = and the function 
(A + X]j=i w i( T ) OA?) ( r '0 vanishes when |r| < |£|. Morever, this shows 
that the Fourier transform of the vector potential A(r, £) is perpendicular to 
the (n + l)-dimensional vector (1,w(t, £)) as 

n 

(4) + X>i( T >*Mj)'W) = (1,o;(t,0)-^(t,0- (42) 

Equation (j4ip has infinitely many solutions and as a matter of fact they can 
be parametrized by S n ~ 2 . On the other hand, equation (j42p tells us that 
^4 = (Ao, . . . , A n ) is orthogonal to all elements of E = {(l,u(r, £)) : r + 
u;(r, £) ■ £ = 0}. It is not hard to prove (see Appendix A) that the orthogonal 
complement E x is one dimensional and since (r, £) is perpendicular to any 
vector of the form (l,u;(r, £)), this complement has to agree with the line 
{c(r,0 : ceR}. 

Since the previous argument works for an arbitrary r and since the set {£ : 
|t| < is an open subset in IR n , we see that vA(r, £) = (v4 (r, £), . . . , A n (r, £)) 
is proportional to the vector (r, £ ) in the complement of the light cone. In 
other words, we can find a function $ such that 

(A (r, 0, • • • , Mr, 0) = i$(r, (r, f) (43) 

whenever |r| < |£|. Since for any j the function Aj decays exponentially in 
t and is compactly supported in x then its Fourier transform A, is analytic 
in the strip |Im r| < c. 

On the other hand, equation (T4"3"|) gives 

= l<j<n, 
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which tell us that $ is analytic in the set {(r, £) : |Im r| < c, (r, £) 7^ 
(0,0)}. Hartog's theorem (see [TT]) tells us that the concepts of removable 
singularities and isolated singularities agree in functions of several complex 
variables and we conclude that $ is analytic in the strip |Im r| < c. Moreover 
if we let ip be the inverse Fourier transform of $, then tp and all of its 
derivatives are exponentially decaying in t and we only need to make sure 
that it has the right support properties. 

Because of the assumptions on the support of the functions Aj we have, 
by the Paley- Wiener theorem 

, T , C^exp^Im^J) 



for 1 < j < n. Then for \t^'\ > 1 



Cfexp^In^) 

- (l+iei)* ' u - J - n ' 



for some CV > 0. Since the function h(r, £) = $(r,£)(l + |£|)^exp(-i2|Im£|) 
is continuous when < i ; it is also bounded, hence \h(r, £)| < Cn for 
some positive Cn and the estimate 



C exp(#|Im e|) 



l^0| <^=^ (44) 



holds for any ( Gl" Making use once again of the Paley- Wiener theorem 
we conclude that the inverse fourier transform of $(r, £) is supported in the 
set {x : \x\ < R}. □ 

Before going ahead to prove the correspoding equality of the scalar poten- 
tials, we will pause for a second to relax the conditions imposed on the vector 
potential. Let us start by replacing exponentially decaying by Schwartz func- 
tions. 

Theorem 6.2. Suppose that A(t, x) = (A (t, x), . . . , A n (t, x)) with A E C°° 
in x and t is such that for any M > and non-negative integers a, (3 there 
exist constants Cm,u,p > suc ^ that (1 + \t\) M \d^d^Aj(t,x)\ < Cm,<x,p f or 
< j < n. If in addition Aj(t,x) = for \x\ > R > and ( fff^j) holds, then 
there exists (p(t,x) G C°°(t,x) such that 
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i) A (t, x) = d t ip(t,x), Aj(t,x) = d x .ip(t,x), 1 < j < n, and 



ii) Suppt^ CRx {\x\ < R}, 



\t\) M \d?dgcp(t,x)\<C' M ^. 



Proof. The proof goes along the same lines as the previous proposition except 
that now in equation ( H3j) we only know that the left hand side is entire in £. 
For To 7^ fixed, Hartog's theorem tells us that $(ro,£) is entire and when 



0, equation (TJ3]) gives $(0, £) 



for 1 < j < n showing that 



$ has no singularities. 

The part of the proof that deals with the support of (p remains unchanged 
and we only need to show that $ is a Schwartz function. If M' > and /3 is 
a non- negative integer, we have for |£| < R 



:i + MH<9f$(r,o| = (i + M) 

= {l + \r\f 
<C{1 + \t\) m 



j=0 



3=0 



< C 



M',[3,W 



where in the last inequality we used the fact that the exponentially decaying 
C°° functions Fourier transform into Schwartz functions. Since $ is itself 
Schwartz, the desired function tp is again the inverse Fourier transform of 
$. □ 

The conditions on the vector potential imposed so far are such that we 
end up working with functions once we compute the Fourier transform of 
equation f )39|) . nevertheless, this transform can be computed under weaker 
assumptions and the following theorem tells us that the final result is still 
valid. 



Theorem 6.3. Suppose that A(t, x) = (A (t, x), . . . , A n (t, x)) with A G C°° 
in x and t is such that for < j < n, \Aj(t, x)\ < C(l + \t\) M with C, M > 
and \t\ > to- If in addition the functions \Aj(t,x)\ are locally integrable in 
R n+1 , satisfy the support condition A j(t, x) = for \x\ > R > 0, and equation 
/ T5P|) holds; then there exists (f(t,x) G C°°(t,x) such that 
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i) A (t, x) = d t ip(t,x), Aj(t,x) = d x .ip(t,x), 1 < j < n, and 

ii) Suppy? CRx {|x| < R}. 

Proof. By the hypothesis on the growth of Aj, 1 < j < n, we can com- 
pute the Fourier transform of equation (|39|) to obtain 5(t + to ■ £)(A + 
Y^ 1 j=i UJ jAj) A ('i~,£) = 0, where Aj(r, £) is an analytic function in £ and a dis- 
tribution in r. In addition, since the wavefront set of 5(t + u • £) and A + 

YTj=i u j ' Aj do not intersect, we can define a restriction of A + YTj=i u j ' Aj 
on the hyperplane r + u ■ £ = (cf. Hormander [12J). Proceeding as before 
we find that when |t| < |£| there are infinitely many solutions of equation 
( I4ip and that they can be parametrized by S n ~ 2 . Moreover, the change 
(r, £) —> (err, a£), a > 0, in (jUJ leads to 

w(ar, «£) 



Mlfl " l a ll r l 

— ■ w(ar,a£) = -— 

which tells us that the solutions co(t, £) of ( )4TT) are homogeneous of degree 
in (r,£). 
Therefore 

n 

A)(r,0 + X) w A( r '0 = (45) 

on the plane r + w • £ = 0. 

Replace (r, £) by (ar, a£) where a > and we then see that for a > 

n 

A (ar, a£) + w^^ar, a£) = 0. (46) 

3=1 

We now let xi a ) be an arbitrary C^°(R) function with support contained in 
the set | a — 1| < e and multiply (f46l) by x( a )- Integration in a leads to 

n 

ao(r,0 + X)w i o i (r,0 = 0, (47) 

i=i 

where r + oj ■ £ = and 

POO 

a i( r >0= / ^jKaOx(a) dot- 
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Notice that Oj(r, £) are no longer distributions and we can put uj = w(r, £) 
in ( HTj) . Arguing as before we find that (cio(t, £), . . . , a n (r, £)) = zfe(r, C)( t j 
for some 6(r, £), or in other words, 

— - — = —i = ■ ■ ■ = — = z6(r, 0- 

Since x( a ) i s arbitrary we get 

A (aT, a£) Ai(a>T,a£) A n (ar, a£) 



where $(ar, a£) is a distribution in ar for all a G (1 — e, 1 + e). Finally when 
a = 1 we get 

A)(t,0 = *t*(t,£), ^i(T.0=i6*(r,0. ••• ,^n(r,0=^n$(r,0 

for |r| < |£|. As before we have that \I/ is entire in £ and f G S' in r 
(since Aj G S" in r). Therefore y? = J^~^ G 5" in £. Moreover the identies 
d t (p = A ,d Xl (p = Ax,...d Xn cp = A n imply that <p(t,x) G C°° in (t,x) and 
that (p = for |x| > _R. □ 

Summarizing, the three previous results prove that the vector potentials 
and are gauge equivalent with gauge g = e ltf> . Next, we show that 
this equivalence implies the equality of the scalar potentials and V^ 2 \ 

By the previous proposition A^ — A^ = V^y?, replacing the pair 
(*4«,1/W) by (A^ 3 \V^) where = A® + V t> ' x cp and V« = we 
find by means of proposition 12.11 that A^ = A^. Next we use our Green's 
formula Q32]) with the pair of potentials (A^, V^) and (.4 (3) , 1/ (3) ) to obtain 

= <(^ 3 > - VV)u,v) [TuT2]xn = £ jT (V^ - W 2 >) uvdxdt. 

Making use of the geometric optics representations (I33|) - (l36l) the above inte- 
gral becomes 



0= f 2 [ (V®(t,aO-V (2) (t,x))e< fla ^ 



Ti Jn 



(48) 
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where • • • denotes terms of order O {k^ 1 ). Taking the limit as k — > +00 we 
notice that the equality of the vector potentials imply that R2 — R\ = 0. 
Thus, after a change of variables, we can rewrite f HHj) as 

= / (7 00 V® (X' + s(l, u)) - (X' + s(l, w)) ds) x 2 {X') dS X ; 

since x is arbitrary the inner integral in the expression above vanishes 

POO 

+ s, x + su) — V^(t' + s,x' + su)) ds = (49) 

»/ 00 

which shows that the light ray transform of the potentials agree. A sim- 
ple variation of the previous proof applies and we have = = . 
Therefore the pair of potentials (A^\ V^) and (JS 2 \ V^) are gauge equiv- 
alent (see also [25], [23]). 

7 Stability estimate 

In this section we assume that the components of the vector potentials A^ 
and A^ are real valued, smooth and compactly supported in both t and x. 
Just as we did before, let us write 

A = A (1) -A {2) where A ik) = (A { *\ . . . , A^), k = l,2, 

and let us further assume that the potential A satisfies the divergence con- 
dition 

•a 

divA = d t A (t, x) + J2 d^Ajit, x) = 0. (50) 

Since the potentials are compactly supported we can find real numbers T\ < 
< T 2 and an open set V C R t x R™ such that if we set Q = {T u T 2 ) x ft, 
then 

i) Q CP; and 

ii) Supp(Aj) C T> for j = 0, 1, . . . , n. 

In other words, T> is a set containing Q and the support of all the components 
of the vector potential. Since T> is bounded we can find T 3 bigger than |Ti| 
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and \T 2 \ such that V C (-T, T) x R™, we then choose T > T 3 + diam(ft) and 
set Q T : = (— T, T) x ft. 

When T is selected appropriately we can find solutions 

u (t,x) = x w (*,x)e ifc( *- w -* )+ifll(t>a!;w) (1 + 0(k~ 1 )) 

and 

v (t,x) = x,(t,i)e ifc(i -" x)+!fl2(M;w) (1 + 0(k- 1 )) 

of the backward and forward hyperbolic equations satisfying 

u = u t = on {-T} x ft, 
V = Vt = on {T} x ft, 

where the function Xu is such that 

a ) + Sj=i Ujdx^Xwit, x) = 0. This is, is constant along light rays; 
and 

b) Xoj is supported in a small neighborhood of the ray {(t + s,x + sou) : 

s e R}. 

As before we can make use of the Green's formula developed in previous 
sections to obtain 

n 

(( A i - A 2)(f),9)^T,T)xdn = J Qt : = i( A i u > ~ id ^ v )q T + ( A i(- id zj u ^ v ) QT ) 

3=1 

n 

+ Y^(((Aff-(Aff)u,v) QT -(Vu,v) QT 

3=1 

- (A u, -idtv) Qr - (A (-id t u),v) QT 
-(((A^-(A^)u,v) Qt 

(51) 

where we have set 

/ = u(t, x) | ( _ TjT)xan # = v(t, x) | ( _ TjT)xan - 
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We next regard Ai — A 2 as a map from H 1 — > H° and denoting by 1 1 j 
the operator norm between these spaces, we have by the Cauchy-Schwarz 
inequality 

\Iq t \ = I ((Ai-A 2 )(/),9) K r)xffl| < |||Ai-A 2 |||x 

WfWrn ((— t,t) xan) \ \g\\ l 2 ((-t,t) x an) ■ 

The latter norm can be estimated by 



II 9 ||l2((-t,t)x£K2) = II Xu(t, + 0(k 1 )) ||i2((_ TjT)xa n) 
< 1 1 Xuj(t, x) I \ L *((-T,T)xdn) + Oik' 1 ), 

whereas the middle norm can be estimated by 



/ ||ffi((-x,T)xafi) < C(n, Q) k\\ Xu \\& {( _ TX)xaa) + 



(52) 



C(n,tt)k 



Xu 



\L 2 ((-T,T)xdn) 



0{k 



(53) 



If in addition we assume that | \xw\ |l 2 ((-t,t)x9q) < C we have by (I52"j) and 

m 



I QT \<kC(n,Q) |||A 1 -A 2 ||| + 0(A; 



-In 



(54) 



where C(n,Q) is a constant that depends upon the supremum of all ray 
integrals of A + YTj=i u j A j- O n the other hand, the integral Iq t leads, just 
as before, to 



/T (■ n 
\ ( A o+y2 u} o A j)( t , x )xi(t, 
-TJn v =1 7 



x) x 



cl.s 



dx dt + 



where "• • ■ " represents terms that, when divided by k, go to zero as k — > +oo. 

Dividing the above expression by k and taking the limit as k — > oo we get, 
after using the triangle inequality and performing the change of coordinates 
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(t, x) = a(l, u) + Y' with Y' e Il {ltU]) 

! ! (A + J2^MY> + a(l,u))xl(Y> 



<C(n,fi)|||Ai-A 2 |||. (55) 



If we set 



if!.. (a,+E"=i<^) (y>+s(1,u)) 



da, 



» n 

a(y') := / (A +^a;^ J )(r+a(l,a;))e- 
j'=i 

equation fl55l) can be rewritten as 

/ a(r) x 2 (r)d^ <C(n,n)ll|Ai-A 2 |||. 
^ n (l,-) 

On the other hand the conditions imposed on the support of \ guaran- 
tee that the above estimate holds true for any such function satisfying the 
condition J n \x(Y')\ 2 dS Y > < 1, thus a is a bouded linear functional on 

L 1 (Il(i iLi; )) and the estimate 



/oo n 
(A Q + Y, UJ j A j ){X' + a{l,oj)) 



x 



e i /-«,(>*o+5^=i « rf Af)(*'+«(i, w )) d, da 



<C(n,fi)|||A 1 -A 2 | 



holds. The Fundamental Theorem of Calculus then gives (in the original 
coordinate system) 



exp 



(t+s, x+stu) ds 



-1 



< C(n,0)|||Ai-A 2 |||. (56) 



In the case of uniqueness of the potentials it was very easy to go from 
an expression concerning the above complex exponential to an expression 
involving only the ray transform of the function Aq + YTj=i UJ j J ^j- m this 
case, obtaining such an estimate is slightly harder and we need to assume 
that the following condition holds: 
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Hi) the supremum 



a :- 



Sn P(t,x;u)&QxS n ~ 



/oo n 
(A + UjAj) (t + s, x + su) ds 



satisfies the inequality a < 2n. 
Denoting by f3 the integral /^(A) + YTj=i u jAj)(t + s, x + so;) ds we get 

le^-ll ' ' ; 



sin 



(57) 



Also, condition in) gives < ^ < 7r and the right hand side of (|57|) is 

bounded from below by a positive constant — . We then have the estimate 

64 



/oo n 
(A +^^ UjAj)(t+s, x+su) ds 
.7=1 



which together with (156]) gives 



< Ca 



e i IZo( A 0+T,j=l WjAj){t+3,X+3Ul) ds_ l 



/OO n 
(A + ^2ujAj)(t + s,x + su) ds < C(n,fi)|||Ai - A 2 | 
3=1 



(58) 



We next want to use fl58l) as well as the divergence condition imposed 
on the potentials to obtain an estimate for the potentials Aj, j — 0, . . . ,n, 
following the ideas in Begmatov's paper [2]. 

If we let F denote the ray transform of Aq + Y^=i along light rays, 
we have 



F : IR f x IR" x S* 1 - 1 -> 



r n 

F(t, x;oj) :— / (A + \^UjAj)(t + s,x + su) ds. 

3=1 

and by (1381 

|F(t,ar;w)|<C7(n,n)|||A 1 -A 2 ||| 
for all (t,x) GR(Xi>G S"" 1 . 



(59) 



(60) 
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The Fourier transform of F in the variables X\ , . . . , x n is 
(J' (x _ +€) F(t,.;w))(0= / e"** / (Aj+rw j Aj)(t + s ) i + si<;)dsdx. 

and the change of coordinates x = x + soj, t = t + s, with Jacobian | | = 1 
leads to 

(F ix ^ F{t,-;u)){0 = e -W [ [ e-^e-^-^HAo + J^^M^^dtdx 

ii" Jm. ^ „-_i ' 



where the right hand side of the above equation is the Fourier transform (in 
all variables) of A + X^Li^'A? a ^ ^ ne P orn t ( — w ' This equation can 
be rewritten as 

n 
3=1 

and we realize that since the right hand side is independent of t, so must be 
the left hand side. In particular when t = we have 

n 

(Aq + J>A) • e,0 = (^5)^(0, •;"))(£) =: Gfcu). (61) 

i=i 

Since the potentials Aj are smooth and compactly supported, F(0, •; •) : 
R™ x S 1 "^ 1 — » R is also smooth and compactly supported^, moreover (1601) 
shows that it is uniformly bounded by C(n, 0)|||Ax — A 2 |||, hence 



|G(£;w)|= / e-^F(0,x;w)dx 



< \\F{0,-;-)\\ LX(M n xS ^)Vol{B n {R)) 

< 67(71,0)^ |||Ai-A 2 |||, (62) 

which tells us that G is uniformly bounded in Rg x S"™ -1 . 

We now turn our attention to the Fourier transform of the potentials. We 
want to obtain an estimate for |A/(t, £)| on a conic set whose complement 



1 This is because for |a;| big enough the light rays with direction (l,w) emanating from 
the point (0,x) do not intersect the support of the potentials Aj. 
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contains the light cone {(r, £) : |r| < |£|} and use this estimate as well as 
an analytic continuation argument to obtain bounds for |-A,-(t, £)| in the full 
space R T x 

For (r, £) fixed with |r| < we know by considerations made in the 
previous section that we can find unit vectors u = w(r, £) parametrized by 
rS n ~ 2 (an (n — 2)-dimensional sphere with radious r, ^ < r < 1), such that 
r + w(r, — an d satisfying w(#r, 6>£) = u(t, £) for any 9 > 0, i.e., w(r, £) 
is homogenous of degree in (r, £). 

We consider a maximal one dimensional sphere with radious r contained 
in rS n ~ 2 and choose unit vectors £) forming the vertices 

of a regular polygon with n sides. We then consider the following set of n + 1 
equations 



j'=i 

n 

(^(r,e) + Ew.e)) = o, 



it 



,n 



(63) 

where the last equation is a simple consequence of the divergence condition 
d t A (t,x) + ^2™ =1 d Xj Aj(t,x) = 0. Our goal is to show that this system is 

uniquely solvable for (A , A 1 , . . . , A n ). 

In order to prove this statement it suffices to show that the matrix 



M(r, 



J ' 

(2), 



r,0 
r,0 



.{A 



(n 



'(r,0 

& 



is invertible. Notice that the entries of M(r, £) are homogeneous of degree in 
(r, £) and the inverse, if it exists, will also have entries that are homogeneous 
of degree in (r, £). 
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To see that M(r, £) is indeed invertible we show the homogeneous system 



A (r,0 +J2uf\r,0Mr,0 = k = 1 

J=1 n (64) 



has no non-trivial solution. Once again, the considerations made in the 
previous section guarantee that the only potentials A = (A Q , A\, . . . , A n ) 
satisfying the first n equations are those of the form A (t, £) = $(r, £)r, 
A?'( r '0 = ^{ T iQ£,ji 3 = 1, for some smooth function <£>. The last 

equation in the above system leads to $(r, £,)\/t 2 + |£| 2 = which in turn 
gives $ = and A = 0. 

Since M(r, £) is invertible we can write 

n 

Mt,0 = Ec M (r,OG(^ W (r,0), 1 < k < n, < j < n, 

k=l 

for some Ckj{r, £) homogeneous of degree in (r, £). This homogeneity prop- 
erty as well as the uniform boundedness of G allows us to compute an estimate 
for the Fourier transform of the potentials Aj in the ray {(ar, ct£) : a G M} 



n 



k=l 
n 

<J2hAr,0\\G{a^^ k \ar,aO)\ 

k=l 

n 



<C(n,0)|||A 1 -A 2 |||^|c fc|i (r,e)|, (65) 

k=l 

where in the last line of the previous inequality we used f )62|) . 

At this point it is convenient to recall that our initial goal is to obtain 
a uniform bound for Aj(r, £) in the set {(r, £) : |r| < y}. I n view of ( 16"5"j) 
it suffices to work on the compact set {(r, £) : r 2 + |£| 2 = 1, |r| < y-}. 
To obtain such a bound it is necessary to study the entries c^r, £) of the 
inverse of the matrix M(r, £). It is a well know result that such entries can 
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be described by 

C ^ (T ^ ) = detM(r,0 C ^ (r ' 
where Cj^r, £) is the (j, fc)-cofactor of M(r, £). 

Since on the set {(r, ^ : r 2 + |£| 2 = 1, |r| < M} all the entries of M(r, f ) 
have absolute value less or equal to one, and since Cj } k(r, £) consists of sums 
of products of n such entries, we have 

Mr,0\< 3t' L < " 



detM(r,OI " |detM(T,{)| 

The quantity | det M(r, £)| can be interpreted as the (n + l)-dimensional 
volume generated by the set of vectors {(1, u^(t, £)), . . . , (1, u/ n )(r, £)), (r, £)}, 
however, since independent of the values (r, £) the vectors 
are chosen to be the vertices of a regular polygon, the n- dimensional volume 
V generated by {(1, £))}, is constant. We then have 

| det M(t, £)| = V x P(r, £) where P(r, £) is the projection of (r, £) into the 
linear subspace generated by {(1, lo^(t, £)), . . . , (1, w^r, £))}. This projec- 
tion is given by C sin where is the angle between (r, £) and said subspace. 
Since the vectors (1, u/ fc )(r, £)), k = 1, . . . , n, are located in the boundary of 
the light cone (i.e., the set {(r, £) : |r| = |£|}), this angle is bounded below 
by f. Therefore on the set {(r,f) : r 2 + |£| 2 = 1, |r| < f } the value 
| det M(t, £)| is uniformly bounded from below by Vsin | and 

Ti 

\c k j{T,0\ < 



Vsin | 



These observations combined with fl65|) give the uniform estimate 

|5(r,£)| <C(n,fi)|||Ai-A 2 ||| (66) 

on the set {(r,f) : |r| < |i }. 

Our next step is to obtain an upper bound for Aj(r, £), j = 0, . . . , n, on 
the complement of {(r, £) : |r| < y}. To obtain such estimate we first fix r 
and compute upper bounds for all lines that pass through the origin and are 
contained in the hyperplane r = r . 

We will consider the case where the line corresponds to the £„-axis, but 
before doing so, we will need some auxiliary results. 
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Lemma 7.1. Consider the strip 



S = {z 



Zi + iz 2 : z\ E 



R 



z 2 \ < 2|r |7r,r ^ 0} 



and the rays 

Pi = {z : -oo < Z\ < -2\t \,z 2 = 0}, p 2 = {z : 2|r | < Z\ < oo, z 2 = 0} 
in the complex plane C. 

If E = pi U p 2 and G = S \ E is the strip with cuts along the rays p\ and p 2 , 
we have 



where w(z,G, E) is the harmonic measure of E with respect to G. 

This statement is a very well known result about harmonic measures, its 
proof is mostly taken from [2J and it is included here for the purpose of self 
contention. 

Proof. For h > 0, the map 



comformally transforms G into the upper half plane H + = {C = Ci + K2 '■ 
(1 G RX2 > 0}. Under this mapping, the interval I = {z : \z\\ < a,z 2 — 0} 
transforms into the imaginary half axis {( : Ci = 0)C2>0}. The boundary 
of G goees into the real axis and the set E = p 1 Up 2 tranforms into the subset 
of the real axis 



E l = {C : Ci < -e a *l\C, 2 = 0}U{C : |Ci| < U2 = 0}U{C : Ci > e^ /h , C2 = 0}. 



Then by the harmonic principle (see jl]), the values of the harmonic measures 
on E, Ei with respect to the sets G, H + agree, this is 



We also know that the harmonic measure on the right hand side of the 
previous equation can be constructed by means of the Poisson integral for 
the upper half plane 




(67) 




w{z,E,G) 



^(C(2),£i,H+). 




(68) 
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where XEi(t) is the characteristic function of E\. 

Since we are interested in the image of / under the map ((z) and since 
this image is precisely the positive imaginary axis, we may assume without 
loss of generality that £ = C2 > 0. From fl68l) we obtain 

/>/ n P tt+n 2/7T C2(exp(a7r//i) - 1)\ 

P7(C(g), g x , H + ) = - I - - arctan — 69 

n \2 {(2) + exp(air/h) J 

Choosing h = art and using the inequality 

aictm |(f^)< arctan (izi) (70) 
(C2) 2 + e 2e2 

we obtain 

2 /tt fe 



, arctan , — < < 1 



and we conclude that 



^ < vj(z,E,G) < 1. 



□ 



Based on this result we want to 'embed' the £ n -axis into said strip and 
use the bounds on the harmonic measure. To do so we realize that since the 
potentials Aj, j — 0, . . . , n, are compactly supported, the functions A,-(ro, £) 
admits an analytic extension in £ n into the complex plane. If we let 

II = {v = (ui, zaj) : v\ e M, M < 2|t |tt, r 7^ 0}, 
q l = { u = : -00 < z/i < — 2|t- (, i/ 2 = 0}, 

92 = {v = {v u v 2 ) : 2|r | < v x < 00, v 2 = 0} 

and restrict ourselves to the £ n -axis (i.e., £1 = • • • = £ n _i = 0), the estimate 
flBTj) leads to 

2 

- < w{v,E u Gi) < 1, 

where Ei = qi U g 2 an d Gi = II \ i^. 

Denoting by Vj(v) = A,(2t , 0, . . . , 0, u), the restriction of Aj to the £„- 
axis, we have by the two-constant theorem (see [IE] Theorem 9.4.5) 

\vj(v)\ < m] Mf (71) 
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where rrij and Mj are the respective upper bounds of the modulus of v(v) 
on the rays qi and q 2 and on the lines {(^1, v^) : v\ G R, u 2 = — 2|r |7r} and 
{(vi, z/ 2 ) : z^i G R, z/ 2 = 2|t |7t}. 

At this point it is worth to point out that the rays qi and q 2 are contained 
in the set {(r, £) : |r| < -y} anc ^ that we have already computed an estimate 
for in that region (equation IB6T) . To compute Mj we resort to the 

identities 

v j (u) = Cfan) [ e-^ 1+W2 ^^(2r ,0,...,0,x n )da; n 
Jr 

with Wj the Fourier transform of Aj in all variables except x n . Next, we 
realize that these functions are compactly supported in x n and the above 
integrand is nonzero only on a finite subset of the real numbers. Hence 

/a(fl) 
-a(Q) 

where a(Q) is a positive number bigger than diam(fi). Since v = v\ + 2^2 is 
restricted to the strip II we have 

e 2|r |7T 

2|r |7r 

and in particular when v is a real number satisfying — 2 1 to I < v < 2|to| we 
have by ( 171]) 

2 I TQ I 7T 2l 

e 3 777, ? 

< C(fi,n) r^- 

2 1 To 1 3 

All this arguments can be carried out to the case where a line is contained 
in the hyperplane r = r Q , passes through the origin but is not parallel to any 
of the axes. Finally using (166]) we obtain an estimate for the Fourier transform 
of the potentials in the set {(t, £) : |t| > — }, namely 

\Mt,0\ < C(Q,n) - 3 111 \7 2 '" 3 . (72) 

From estimates ( )66|) and ( 1721) we can establish the desired estimate for the 
vector potentials. The general idea is to use the inequality < 11/ IU 1 

and partition R T x RJ? in an appropriate way. 
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From the Fourier inversion formula we have 

A 3 (t, x) = C(tt) f [ S^-VAjir, drd£ (73) 

J </R T xR™ 

and by taking absolute values we have for any pi > 

\Mt,x)\ < C(tt) // |5(r,0|drde 

< C(tt) // |5(r,0|drde 
JJB( P1 ) 



+ C(tt) // |^(r,0|drde 
J JB( Pl y 

= h + h, 

where B(pi) denotes the (n + l)-dimensional ball B(pi) = {(r, £) : |r| 2 + 
l£| 2 <Pi}-' 

To obtain a bound for J 2 we recall that the potentials A,-, j = 1, . . . ,n, 
are C°° in t and x. Hence for any (5 > and any p x > if |r| 2 + |£| 2 < p\ 



(M 2 + I£l 2 ) f 

If /3 > n + 2 the integral J 2 converges. Moreover, the estimate 

,Mr.O\^<^<^ (74) 

B(pi) c Pi Pi 

holds. 

To estimate I\ we break up the ball -B(pi) into two smaller pieces 



Ci = fl(pi)n{(T,0 : M<^} and C 2 = n {(r, : M > ^=}. 

Then 



/i< // |^(r,0|drd£+ // |^(r,0|drde, 
and since Ci is a compact subset of B(pi) we have 

/i<Cp" +1 + // |Aj(r,0|drde 
35 



The advantage of this decomposition is that C 2 is contained in the set {(r, £) : 



t\ > -y} and that on this set \t\ is bounded below by -^=. Thus by ((72 



Pi 



75' 



(75) 



where we have set p 



Pi 



Equations (E3])-(J75]) lead to 



\Aj(t,x)\ < C(tt,n) 
< C(Q,n) 



+ p n+1 +p n+ t e ^|||A 1 - A 2 



2 2 P 7T 



-1 
-p 



. 2 2p;r 



+ p n+ 3 e~3~ | | | Al - Aj 



(76) 



Now the idea is to choose p small enough so that the two terms in the the 
right hand side of (J76l) are comparable. In other words we want p to satisfy 
the identity 



2 2pn 



c 

— = p" ' 3 e s 
P 



A, 



for some constant C. By taking logarithms on both sides of the previous 
equation we obtain the equivalent identity 

C 



2 log 



A, 



(3n + 5) logp + 27rp, 



(77) 



and since the right hand side of ( 1771) is one to one when p > we know that 
it admits a unique solution. 

On the other hand, the inequality log p < p for positive p as well as (J77J) 
lead to 

C 

2 log — j-fjy < (3n + 5 + 2vr)p, 



or 



i - A 2 



1 3n + 5 + 2vr 
p- 2 



log- 



C 



|Ai-A 5 



and equation (jTbj) becomes 



log 



C 



i - A 2 | 

Summarizing, we have proved the following 



36 



Theorem 7.1. Suppose that the vector potentials A® = (Aq\ . . . , An), I = 
1,2, are real valued, compactly supported and C°° in t and x. Let A = 
(Aq, Ax, . . . , A n ) where Aj = Aj — Aj and suppose that the divergence 
condition 

n 

divA = d t A (t, x) + J2 dxAj(t, x) = 

3=1 

holds. If A; represents the Dirichlet to Neumann operator associated to the 
hyperbolic problem (l)-(4), then the stability estimate 

-i 



max 

0<j<n 



Af\t,x)-Af\t,x) 



L°°(R t xR5) 



log 



c 



|Ai- A 2 



(78) 



holds for Ai, A 2 satisfying |||Ai — A 2 ||| << 1. 



Proof. The hypothesis of the theorem guarantee that conditions i) and ii) 
hold and the only condition that is not automatically satisfied from the hy- 
pothesis is condition Hi). However a simple rescaling of the vector potentials 
Aj — > A'j = \Aj, with a the supremum of all ray integrals of the potentials 
as well as a similar rescaling of the coordinate axis show that the estimate 
(I78p holds for the potentials A'j. In turn, this implies a similar estimate for 
the original potentials Aj. □ 



8 Presence of obstacles 

One variation of the above problem consists in the introduction of convex 
obstacles inside the domain Q. That is, let fi^, 1 < A; < M, be simply 
connected bounded domains in M n , n > 3 and let D — Q \ Ujj^f^. If we 
consider again the equation 

n 

(-td t + A (t,x)) 2 u-^2(-id Xj +Aj(t,x)) 2 u + V(t,x)u = inlxfi, 

3=1 

u(t,x) = for £ << 

u(t, x) = f (t, x) onlx dfl, 

with the additional condition 

«(^)L^ = °> ±<k<M. (79) 
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Then we can prove as in the case of no obstacles that the vector valued 
potential A(t, x) satisfies 



VA(t,x;u) = 




for any ray 

7 = {(t,ar) + s(l,o;)|seM} (81) 

not intersecting the obstacles R x Qk, 1 < k < M. We will show that under 
some conditions the vector and scalar potentials can be recovered outside 
these obstacles. 

As a warm up let us consider the case where n — 3, there is only one 
obstacle and the integrals over light rays are zero for a smooth scalar function 
/ = f(t,x) satisfying the support condition f(t,x) = when \x\ > R, and 
the growth condition \f(t,x)\ < C(l + 1^1)^ for some integer N. 

Under these settings, for x an arbitrary point in D, we can find a two 
dimensional plane P in M 3 containing x such that P misses Oi, this is P fl 
Qi = 0. We then realize that the set of rays 7 of the form ( IHTj) that are 
contained in the three dimensional space RjxP and pass through x, can be 
parametrized by the one dimensional sphere S 1 . This in turn, allows us to 
use our previous considerations for a scalar function in two dimensions (c.f. 
Ramm-J. Sjdstrand, [23]) to conclude that / vanishes in M.] x P. Since P 
could be any two dimensional plane not intersecting Qi and x was selected 
to be an arbitrary point in D = Q \ Q% we conclude that / = on D x Ij. 

Unfortunately for the vector potential things are going to be slightly 
harder and we will have to make some geometric considerations in order to 
obtain an equivalent result. 

We shall impose the following restriction on the problem: 

(Gl) The obstacles are convex and 

— when n > 4, for each x G D there exists a two dimensional plane 
P passing through x such that P does not intersect any of the 
obstacles. 

— when n = 3 there is only one obstacle. 

Theorem 8.1. // condition (Gl) holds and the Dirichlet to Neumann op- 
erators are equal (Ai = A 2 ), then the potentials A^ 1 ' and A^ 1 are gauge 
equivalent, i.e., there exists a smooth function (p such that A^ — A^ = dp. 
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Proof. Let us proceed first by assuming that n = 3. 

Regarding the vector potential A(t, 1-form, the goal will be to 

prove that the 2- form dA vanishes outside Q\. 

For x not in Q\, using condition (Gl) we can find a 2-dimensional plane 
Pclj not intersecting the obstacle and we can choose three linearly inde- 
pendent unit vectors r)j, j = 1,2,3 close to P such that any plane Pj P , 1 < 
j,p < 3 passing through x + rjj and x + r] p does not meet Qi- Next, introduc- 
ing coordinates x' = (x^, x 2 , x 2 ) in M.^ by the formula x — x = x' l r]i + ■ ■ ■ + x' 3 q 3 
and denoting by A' and B' the 1-form and 2-form mentioned above expressed 
in the new coordinate system, we have 

A' = A' dx' Q + A^dxi + A' 2 dx' 2 + A' 3 dx' 3 (82) 
B'(t,x') = b'j p (t,x')dx'j Adx' p (83) 



0<p<i<3 



where x' = t and 



dA 3 dA' p 
bjp= ~dx^~~dxJ' (84) 
The restriction of B' to the 3-dimensional space x Pj p is given by 

B' \ R i xPjp = b' Qj dx' Q A dx'j + b' 0p dx' Q A dx' p + b' jp dx'j A dx' p (85) 
and since there are no obstacles in M.] x Pj p theorem (16. 3p gives 

A '\rIxp jp = V^ p (t,4,x;) (86) 

where V is the gradient in the (t,x'j,x p ) and ipj p is a smooth function 
compactly supported ). Clearly then equation (184p gives 

As the above discussion holds for all three dimensional spaces parallel 
and close to Pj P we see that B' and thus B = dA vanishes for x near x and 
for all values of t. Being that x is an arbitrary point not in the obstacle, we 
get that d*4 vanishes outside Qi and thus, since we are working in a simply 
connected domain, the vector potential is the gradient of a smooth function. 

When there are two or more obstacles and a; is a point not lying in any 
of them, we resort to the geometric condition (Gl) to find a two dimensional 
plane Pq such that Pq intersects no obstacles and choose n linearly indepen- 
dent unit vectors rjj, . . . ,rj n close to Pq in such a way that r/i, r\ 2 6 P and 
the planes Pj P , 1 < j,p < n, do not intersect any obstacle. 
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As before, we introduce new coordinates x' = (x[, . . . , x' n ) given by x—x = 
x 'iVi + • • • + x'nVn, and express A and B in terms of these new coordidates. 
If we now consider the restrictions of B to the spaces M.j x Pj p , 1 < j,p < n, 
we realize thet we are back into the previous case. Making use of the fact 
that ^'| R i xP . = V<^jp(t, x'j, x' p ) we obtain via equation (184]) that 



B 'k«P JP = ^ (87) 

which as before leads to B — and hence A(t, x) = Vt, x ^{t,x) for some 
smooth □ 

When condition (Gl) fails to hold, we can still recover some information 
regarding the difference of the vector potentials provided that d^4 = 0. As 
we did before, let us regard vector potentials as 1-forms and let us consider 
the case when we have 2 obstacles inside the domain Q. 

This time it might be the case that the domain is not simply connected 
and that for some close path 7 the integral J^A-dx is not zero (here we set 
x = (t,x)), however, the condition dA = guarantees that it only depends 
on the homotopy class of 7. 

With this consideration in mind we want to be able to "span" every 
possible homotopy class representative 7 of the domain Q by using light 
rays. In other words, we would like to impose on our domain the geometric 
condition (G2): every homotopy representative can be continuosly contour- 
deformed into light rays. 

If this geometric condition is met, we can write for 71 an arbitrary simple 
closed path in Q surrounding the first of the obstacles 



A(x) ■ dx = / A(x) ■ dx + / A(x) ■ dx H h / A(x) ■ dx, 

where the set of light rays £1,1%, . . . ,£ r , surround the first of the obstacles 
and are such that they do not intersect the second obstacle. Then by the 
previous arguments regarding the construction of geometric optics solutions 
and Green's formula we have if A^ and A^ correspond to equal Dirichlet 
to Neumann operators 



(A {1 \x)-A (2) (x)) ■ dx = 2n mi 



11 
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Proceeding in a similar fashion, we have for the second obstacle and a contour 
72 surrounding it 

/ (A (1 \x) - A (2) (x)) ■ dx = 2nm 2 . 

J 72 

Then, for an arbitrary closed contour 7 we have (after a contour deformation 
argument) that 

/ (A {1) (x) - A (2 \x)) ■ dx = 2nc 1 m 1 + 2Tcc 2 m 2 . 

x'7= c l7l+ c 272 

where c\ and c 2 are two integers. We now let 6 fx), be the function that 
computes the angle between the projection of x into the hyperplane t — 
and the vector (0, 0, ... , 0, 1), and for j — 1, 2 we set ®j(x) = Q(x — pj), 
where pj is any point inside the obstacle j. Then the functions Gj, compute 
the 'angle that a vector makes inside the obstacle j' and we have 

(A {1) (x) - A {2 \x) - miQ^x) - m 2 e 2 (x)) • dx = 

for any closed contour 7. Therefore 

A [1 \x) - A {2 \x) - mxQi{x) - m 2 e 2 (x) = d^{x) (88) 

for some function (p. 

We can certainly say more, if 7(^0; x) is any curve joining the points 
x = (t,x) and x = (t ,x ), and we let 

Cq(x) = exp I — i / (m!6i(x) — m 2 9 2 (x)) • da; 

V J-y(x ;x) 

then 

-——d s C {x) = miQi(x) -m 2 6 2 (x) 
C (x) 

and equation (!88|) can be reewritten as 

A w (x) - A^{x) - -^—d s C {x) = d s <p(x) 
C {x) 

or 

A^(x)~A^(x) = -^-d s C(x) 
C{x) 
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where C(x) = C (x) exp(i(p(x)). 

To conclude this section let us consider the case when we have any number 
of obstacles, d.4. = 0, A = A^ — A^, f2 is multi-connected and condition 
(G2) holds. 

Theorem 8.2. If dA^ — dA^ = 0, condition (G2) is satisfied and the 
Dirichlet to Neumann operators are equal (h\ = A.2). Then the potentials A^ 
and A^ are gauge equivalent, i.e., there exists a smooth function C — C(x) 
such that \C\ = 1 for x = (t, x) e R x dtt and A {1) - A i2) = ^d s C(x). 

Proof. Let G({-T,T) x ft) denote the gauge group corresponding to the set 
(— T, T) x Q and let £i, . . . ,£ p be a basis for the homology group, this is, 
7 = n\i\ + • • • + n p £ p for any closed contour 7. Since d.4. = the integrals 
/ A(x) ■ dx depend only on the homotopy class of 7. By condition (G2) 
the basis of the homology group can be spanned by light rays that do not 
intersect any of the obstacles, hence if A^ and A^ correspond to gauge 
equivalent Dirichlet-to-Neumann operators, we have as before that 

Therefore for any closed curve 7 

(.4 {1) (x) - A {2) (x)) ■ dx = 2vrg, q E Z, 

which in turn implies the existance of a function C(t,x) G T, T) x fi) 
such that 

A^(x)-A^(x) = -^-d,C(x). 

□ 



42 



Appendix A: Auxiliary results 

Orthogonal complement of a set in M. m 
Lemma: The orthogonal complement of the set 

E = {(l,u) : ueS( n - 1 \T + u.Z = 0,\T\<\Z\} 

is a one dimensional subspace ofM. n+1 . 

Let us start off with some linear algebra facts. 
Let m > 2. If A C B C R m with 

r 

AC5C Span(A) = | ^ « p a p : a p E R, a p E A, r E n| , 

p=i 

then 

A x = £ x = Span(A) ± . 

Indeed, since orthogonal complements reverse inclusions we have Span (A) - 1 C 
B ± C A- 1 . Also, if rr • a = for all a E A, then x ■ Yl P =o a p a p = J2 P =o ®p(% • 
a p ) = 0, which shows that A 1 C Span(A)- 1 . Therefore A- 1 = B 1 - = 
Span(A) ± . 

Denoting by CH(A) the convex hull of A and by C(A) the cone spanned 
by A we have 

r r 

CH(A) = | ^2 a P a p '■ ^2® p = 1,0 < a p < l,a p E A,r EN^ 
p=i p=i 

C{A) = {ta : t E IR + , a E A}. 

Clearly A C GH(A) C C(CH(A)) and since both sets contain particular 
linear combinations of elements of A we also have Span(c(CH(A))) = 
Span(;4), hence, 

Span(c(CH(A))) = Span(A) ± . 

We want to apply these remarks to the set E but before doing so let us 
recall that for |t| < |£| the vectors u satisfaying \oj\ — 1, t + uj-£, — 0, can be 
parametraized by S n ~ 2 . Since rotations are non-singular transformations we 
can compute instead the dimension of the orthogonal complement of the set 

E = {(1, wi, . . . , w„_i, a) : w 2 H h w 2 ^ = 1 - a 2 }, 
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where < a < 1 is a fixed number. Taking into account our previous 
observations we then have 

E L = Span(c(CH( J B))) ± 

= Span(c({(l,wi, . . . ,£*; n _i, a) : w 2 H hw^ < 1 - a 2 })) , 

= Span({(t,t0 1 ,...,te n _ 1 ,ta) : *A, . . . , n -i e M,# 2 + ••• + ^ < 1 - a 2 })" 

and we can see that -E -1 is a one dimensional subspace of IR ra+1 since clearly 
Span(C(CH(£))) is n-dimensional. 
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